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I — I Abstract 
O 

According to the standard ACDM model, the accelerated expansion of the Universe will go on 

forever. Motivated by recent observational results, we explore the possibility of a finite phase of 

Q acceleration which asymptotically approaches another period of decelerated expansion. Extending 

;_( 

an earlier study on a corresponding homogeneous and isotropic dynamics, in which interactions 

c3 

between dark matter and dark energy are crucial, the present paper also investigates the dynamics 
of the matter perturbations both on the Newtonian and General Relativistic (GR) levels and quan- 
tifies the potential relevance of perturbations of the dark-energy component. In the background, 

cn 

lO the model is tested against the Supernova type la (SNIa) data of the Constitution set and on the 

perturbative level against growth rate data, among them those of the WiggleZ survey, and the data 
I of the 2dFGRS project. Our results indicate that a transient phase of accelerated expansion is not 

excluded by current observations. 
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I. INTRODUCTION 



It is generally believed that our presently observable Universe is dynamically dominated 
by a dark sector which is composed of a dark-energy component with a large negative pres- 
sure and pressureless dark matter. The physical nature of both these components remains a 
mystery, notwithstanding the intense research activity in the field since the discovery of the 
accelerated expansion of the Universe in [1]. Direct and indirect support for this result has 
been accumulating over the past years. This comprises further results from advanced data 
sets for supernovas of type la (SNIa) as well as results from studies of the large scale struc- 
ture [2], cosmic microwave background [3j, the integrated Sachs- Wolfe effect baryonic 
acoustic oscillations [5] and gravitational lensing [6]. By now, there exists a standard model, 
the ACDM model, which, grosso modo, is compatible with the cosmological data. (Notice, 
however, that there is an ongoing discussions, see, e.g., [?], on apparent shortcomings of this 
model). Nevertheless, because of the cosmological constant problem in its different facets, 
including the coincidence problem, a still growing number of competing models has been 
developed over the last years, most of them "dynamizing" the cosmological constant or even 
generalizing Einstein's theory. Observations force these models to have a dynamics that is 
very similar to that of the ACDM model around the present epoch. Moreover, the past 
evolution is restricted by the necessity of a matter-dominated epoch to guarantee cosmic 
structure formation. The future cosmological evolution within alternative models, however, 
may be very different from a de Sitter phase, which is the final fate of a ACDM universe and 
also of several other approaches like Chaplygin-gas scenarios. Phantom-type dark energy 
with a constant equation-of-state parameter, e.g., will end in a big-rip singularity after a fi- 
nite time [8l|9]. More recently, a still different scenario, called "little rip", was proposed [TO] . 
Already at the beginning of the past decade several authors discussed the possibility that the 
currently observed accelerated expansion might be a transient phenomenon, i.e. that there 
might occur a transition back to decelerated expansion [TTHI3]- Some recent observations 
seem to back up this idea. Evidence was found for a slowing-down of the expansion rate 
of the Universe, equivalent to an increase of the deceleration parameter q{z) for decreasing 
redshifts z close to the present epoch z = [TH4T7] . This could indicate a scenario, according 
to which the observed accelerated expansion of the Universe is a transient phenomenon, im- 
plying a transition back to a decelerated expansion either for the future evolution with z ^ 
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or even around the current epoch z ^ 0. Scenarios of transient acceleration were previously 
discussed in |18j and [IH]. The model on which the present paper relies was developed in 
PU] . It describes transient cosmological acceleration as the consequence of an interaction 
between dark matter and dark energy. Such a dynamics cannot be obtained if the interaction 
represents a small correction to the standard ACDM model. For models of such type the 
long-time cosmological dynamics will always be determined by the cosmological term and 
result in accelerated expansion. To achieve transient accelerated expansion, a twofold role of 
the interaction is necessary. At first, it has to cancel the "bare" cosmological constant and 
at second it has to generate a phase of accelerated expansion by itself. Acceleration has to 
be an interaction phenomenon. As it was shown in [20], these requirements can be fulfilled 
by interaction terms that combine powers and exponentials of the cosmic scale factor. Even 
though this specific choice was made for mathematical convenience, we expect that the men- 
tioned two features will be crucial for a broader class of models. While the study in [20] was 
restricted to the homogeneous and isotropic background, the present paper investigates the 
corresponding perturbation dynamics as well. In particular, we calculate the growth rate of 
the matter perturbations and the matter power spectrum. We compare our results with the 
growth-rate data collected in [21] as well as with those of the WiggleZ survey [22] and with 
the data from the 2dFGRS program [23]. Emphasis is also put on the potential relevance 
of perturbations of the dark-energy component which are neglected by many studies of the 
matter perturbation behavior (for exceptions see, e.g. [2H - EB] ). In concordance with parallel 
investigations for other models [27], we find that they are small indeed on scales that are 
relevant for structure formation. However, there are indication that their role is increasing 
with increasing scale. The background dynamics is reconsidered on the basis of the SNIa 
observations by the Constitution set |28j . 

The paper is organized as follows. In section |Tl] we reanalyze the basic features and 
the homogeneous and isotropic background dynamics of the transient acceleration model. 



Section |III| is devoted to a Newtonian treatment of the perturbation dynamics. A fully 
relativistic, gauge-invariant investigation and a calculation of the matter power spectrum 



are the subjects of section IV In section M we summarize and discuss our results. 
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II. THE TRANSIENT ACCELERATION MODEL 



We assume the cosmic substratum to be dynamically dominated by a mixture of a pres- 
sureless matter fluid and a dark-energy component. The fleld equations for a spatially flat, 
homogeneous and isotropic two-component universe of this type are the Friedmann equation 

3ff2 = 87rG (p„ + p,) (1) 

and 

H = -AttG (Pm + Pcc+Px) . (2) 

Here, pm is the energy density of pressureless dark matter and px is the density of the dark- 
energy component with a pressure Px- The Hubble rate H is given by if = ^, where a is the 
scale factor of the Robertson- Walker metric and a dot denotes the derivative with respect 
to the cosmic time. We assume that both the dark components do not conserve separately 
but interact with each other in such a manner that the balance equations take the form 

Pm + SHpm = Q (3) 

and 

Px + 3H{1 + w)px = -Q , (4) 

where w = ^ is the equation-of-state parameter of the dark energy. The sum of (|3| and 
(|4]) results in the total energy conservation equation p + 3H (p + p) = 0, where the total 
pressure equals the dark energy pressure, p = px- 

It is convenient to write the energy density of the matter fluid as 

Pm = PmoOT^ f (a) , (5) 

where we have chosen = 1 for the present value of the scale factor. The quantity p^o is 
a constant and /(a) is an arbitrary time-dependent function. With / (a) = 1 + g {a), this 
structure implies that 

Q = Pm^ = pmoa~^ f = Pmoj^aa-^ (6) 

and 

pm = pmo (1 + g) a'^ ■ (7) 
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The present values and are related by 



mo 



(8) 



where go = (?(!). The quantity pmo is the value of pm at a = 1 in the presence of the 
interaction, pmo is the value of pm at a = 1 for vanishing interaction. The interaction 
re-normalizes the present value of pm- 

In |20] it has been shown that an analytically solvable transient acceleration scenario can 
be based on an equation-of-state parameter w = —1 with an interaction, characterized by 



g{a) = ca exp(— a /o" ) . 



(9) 



In the following, we briefly recall the basic features of this approach. We start by integrating 
equation ^ with ^ and ^ which yields 



Pm = Pmoa ^ + Ka ^ [a^exp(-a^/(T^) - exp(-l/cr 
where K = cpmo, while from Q, ^ and (|9| it follows that 

Px = P"JJ - K exp {-a'la^) 



(10) 



and 



P"JJ = Pxo - exp(-l/(T2) 



(11) 



(12) 



In the interaction-free limit i^' — t- we have consistently px — )■ Pxo = const. The quantity 
p^^-^ can be seen as an effective cosmological constant which is re-normalized compared with 
the "bare" value, corresponding to p^o, due to the presence of an interaction. The ratio | 
becomes 



a 
a 



Vtmo - is:exp(-l/(T2) 



2Q.l{f + 3irexp(-aV(T2) [a^ - a^] 



with K = and il^^J = The present value of the deceleration parameter is 

a , 1 



(13) 



{i + 3wnxo} 



(14) 



where fl 



XO 



To have a viable cosmological model of transient acceleration, formula (13) should admit 



a transition from ^ < to ^ > before the present time, i.e., for a < 1. If, moreover, the 
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accelerated expansion is a transient phenomenon, there should be a change back from ^ > 
to ^ < at some time, which may be close to the present epoch or at a future period a > 1. 



In the expression (13) the a terms on the right hand side dominate for small values of a, 



i.e., there is decelerated expansion for a ^ 1 provided the condition 

Qmo > Kexp{-l/(r^) 



(15) 



is satisfied. This condition puts an upper limit on the admissible interaction strength. In 
the non- interacting limit it just expresses the positivity of the matter energy density. Let's 
consider now the case a ^ 1. The dominating contribution in the braces on the right-hand 



side of (13) are then given by the constant term —Ip^JJ. As long as p^JJ > 0, however small 
it may be, we will have f > for a ^ 1, i.e., there is no transition back to decelerated 
expansion. This holds, in particular, in the non-interacting limit which reproduces the 
ACDM model. Then Q'^^J reduces to Qxo^ equivalent to ^Iao- Fo^' o ^ 1 this term will 
always dominate the dynamics. An obvious way to obtain decelerated expansion for a ^ 1 



is to put p 



eff 



in (13). This corresponds to a vanishing total cosmological constant. 



In other words, part of the interaction cancels the "bare" cosmological constant, described 
by pxQ. Under this condition it is exclusively the remaining part of the interaction which 
potentially can trigger a period of accelerated expansion. In such a case one obtains from 



(11) that 



n 



xo 



K exp(-l/(T2) 



-a 



1 



(16) 



Then the energy densities of the dark components are (10) for pm and (11) for p^ with 
PxQ = 0, I.e., ^ 

p. = ^a^irexp(-aV^^)(l-^^) • (17) 

Notice that for > a positive value of p. requires a < ^fa. For any a > a the entire 



quantity (17) is exponentially suppressed, the amount of p^ tends to zero (possibly through 



an intermediate period with p^ < 0). With p.^ from (10) and p^ from (17), the Hubble rate 
for this model is 



772 



h -cr A exp(-a a ) . 

a-^ 2 



(18) 



Both for a ^ 1 and for a ^ 1 one has 



h 0.65 




0.030 
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FIG. 1: Two-dimensional probability contours {la, 2a and 3a), based on the Constitution data, 
for all combinations of the parameters h, K and a. 



/Vmin 


a 


K 


h 


465.5 


tr r,o+0.05 

^•^"-*-0.05 


n ni O+0.0004 

U.Ui8_o.0004 


'-'•"^-0.003 



TABLE I: Best-fit values, based on the Constitution data, for the parameters h, K and a. 



The acceleration equation becomes 
1 



a 
a 



1 - li^cr^expf-l/o-^) - , 
? — - + 3irexp(-aV(T2) [a^ - a^] 



To have decelerated expansion for a <^ 1, 



Ka'^ exp(-l/(T2) < 



(20) 



(21) 



has to be required. This condition is similar to (15). The zeros of (20) determine the values 
aq of a at which transitions between decelerated and accelerated expansion (or the reverse) 
occur, namely 
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-^K exp(-l/(T2) + exp(-aj/a2) [a^ - aj] = 1 . 



The condition to have acceleration at the present epoch with a = 1 is 



(22) 



am 



|o > ^ K exp(-l/(T2) 



> 



(23) 



If the inequality (23) holds, we may have present acceleration under the condition pl^/ = 0, 



i.e., a vanishing total cosmological constant. Obviously, the normalized interaction strength 
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K has to be larger than a threshold value to realize this configuration. The condition (23) 



is consistent with (14) if the latter is combined with (16). On the other hand, we have the 



upper limit (|21|). This means, there exists a range for admissible values of the interaction 

2 e^l'^' 



strength, determined by 



9.^-1 



<K < 



2e^ 
30-2 



(24) 



The parameters K and cx^ enter the present ratio of the energy densities for which we find 

M = ^exp(-l/a2)(|a2_i) 

Pmo l-irexp(-l/a2)(|a2-l) ■ ^ > 

The results of a Bayesian statistical analysis, using the 397 SNIa data of the Constitution 
sample [28], are shown in Fig. [l] and in Table [l| For the best-fit values in Table [l] the 
inequalities (24) are satisfied and for the ratio (25) we obtain = 2.275]'^q 285. '^^^ range 
in (24) specifies to 0.0086 < K < 0.0253. (For the integration a range -0.1 < ^ < 0.1 
was used but robustness tests showed that the results are independent of this choice.) The 



numerical value for the left-hand side in (21) becomes 0.475 < | and for the right-hand side 



of (23) we find 0.463 > 



9' 



The relation (20) may be compared with the corresponding expression of the ACDM 
model: 



I ACDM 



2Qf 



(26) 



The interaction term in (20) plays the role of Qa in (26). Now we know that the ACDM 
model provides a fairly good description of the present universe, i.e., for a = 1. This suggest 
positive values of the interaction constant K together with a > 1, which is indeed confirmed 
by our analysis. 



The relation between p^o in (10) and pmo in the definition of K following (10) is 

Pmo 



Pmo 



With the help of the definition Q. 



mo 



1 + c exp (— l/(j2) 

SttGp^q the constant c may be written as 
K 



(27) 



VLrr,o - K exp (-1/0-2) 
The quantity in (|9]) is then given by 

K exp (-a2/o-2) 



fi^o - exp (-1/0-2) 1 
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K exp (— a2/o-2) 
- la'^K exp(-l/(T2) 



(28) 



(29) 



and the interaction term — = t4- becomes 

Pm 1+9 



Q - f5-2f^ja5exp(-aVa') g f a'\ 

P,n " - K [exp (-1/^2) - a5 exp (-aV^^)] " 1 + V ' 

It is obvious that a transfer of energy from dark energy to dark matter, characterized by 
Q > 0, requires < |(T^. While positive values of Q seem to be favored on thermodynamical 
grounds [29j, the observational situation is less clear [30]. Moreover, negative values of both 



the energy density (17) and of Q in (30) for sufficiently large values of the scale factor are 
exponentially suppressed and the total energy remains always positive. According to this 
scenario dark energy is transformed into dark matter in the past, at present and in the 
future untU a < ^ . S. M a ^ . e the deus.ty of the da.. e„e,, beco.es 

negative. However, it decays only until a ^ 8. At this point the direction of the energy 
transfer reverses and for a ^ a it tends to zero exponentially. 
The background interaction term Q may also be written as 

Q = SfiHp^ = fiQp^ (31) 

with 

2 5 — 2 — 

p{a) = -^^. (32) 

^ a2 3 

The consistency of the model can be checked by realizing = (cf.Eq. (|4])) for w = —1) 
with px from (17) and Q from (30) with pm from ([T]) with g from The behavior of the 



deceleration parameter for the best-fit values of Table [T] is shown in Fig. |2} According to our 
model, the Universe is still in accelerated expansion at the present epoch but q{z) will go 
through a minimum in the future and enter a phase of decelerated expansion again. In the 
far-future limit a ^ 1 (not shown in the figure) q will approach q = \ again. 



III. NEWTONIAN PERTURBATION THEORY 



In this section we focus on perturbations with wavelengths much smaller than the Hubble 
radius. Under this condition, the dynamics is well approximated by a Newtonian analysis. 
Afterwards we shall clarify how this approximation fits into a general relativistic scheme. 
Starting point for a Newtonian treatment is the matter energy balance equation + 
(pm,f") a~ Qi where greek indices run over 1, 2 and 3 and is the (non-relativistic) matter 
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1.0 



q(z) 0.0 




FIG. 2: The deceleration parameter of the transient acceleration model as function of the red- 
shift for the best-fit parameters in Table [l] (solid line). The dashed line shows the corresponding 
dependence for the ACDM model. The value 1 = \ corresponds to the Einstein-de Sitter universe. 



velocity. The perturbed energy balance is, in first order, 



pm + + Pmv'l, = Q . 



(33) 



Here, a hat on top of the symbol denotes the first-order perturbation of the corresponding 
quantity. It i 
equivalent to 



quantity. It is convenient to introduce the fractional quantity Sm = Eq- (|33|) is then 



Pr, 



Q — Q^n 



(34) 



The right-hand side of this equation describes the influence of the interaction on the pertur- 
bation dynamics. In the interaction-free limit it reduces to zero. Since the Newtonian model 
does not specify Q we shall assume for simplicity Q = PQ6m. where /3 is a constant. In the 
subsequent section we shall look at this term more carefully. The limit /3 = corresponds 
to an interaction only in the background. For /3 = 1 the interaction does not affect the 
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perturbation dynamics. Eq. (34) specifies to 

+ =-(1 



/3) . 

Pm 



(35) 



Assuming also, that there is separate momentum conservation of both components (the more 
general case of a coupling also via exchange of momentum will be considered in the following 



section), the non-relativistic Euler equation for the matter reads ^ + {v^V 
where is the gravitational potential. From the first-order Euler equation we find 



--Vn 



(36) 



Introducing comoving coordinate by = ag", differentiating (36) with respect to and 



combing the result with (35|) results in 



Q 



Pr. 



Sm + {1- /3) 



Q 

2H— + 

Pm 



Q_ 

Pm 



(37) 



where is the Laplacian with respect to the comoving coordinates. Equation (37) demon- 
strates the influence of the interaction on the perturbation dynamics. Both the coefficients 
of 5m and 5m depend on Q explicitly. Even for /3 = 1, the case in which the interaction is 
not directly felt at the perturbative level, the Hubble rate H is essentially determined by 



the interaction according to (18). The first-order field equation of Newtonian gravity reads 

1 



AttG {5pm + 5px) = AttG {pm5m + Px5x) 



(38) 



where 5x = In many studies of the growth rate of matter perturbations the dark-energy 



perturbations are neglected. This corresponds to assuming 5^; = in (38). However, this is 
strictly justified only for a cosmological constant. In dynamical dark-energy models 5^ is dif- 
ferent from zero and the matter perturbations are coupled to the dark-energy perturbation. 
Neglecting this influence may result in an incorrect interpretation of observational data . 
On the other hand, for specific models the coupling can indeed be shown to be negligible on 
small scales [27]. In order to obtain a closed second-order equation for 5m we shall assume 
here a simple relation of proportionality 5x = ot5m between 5x and 5mi where a is constant 
(cf. [31]). The limit a = corresponds to vanishing dark-energy fluctuations. For any a of 
the order of one, the dark-energy perturbations are relevant for structure formation. Under 
this condition we have 



AttG {pm + apx) 5„ 



a) h a 

P 



5r, 



(39) 
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for the last term in Eq. (37). For the term that multiphes 5m in (37) one finds 

1 



Q_ 

Pm 



Q 

2H — = H 

Pm 



2 5-24 



Pm 



where 



with 



and 



1 - la^K [exp {-l/a^) - exp {-a^/a^)] 



Bio) = 



5-24 



(40) 

(41) 
(42) 
(43) 



With the help of these relations and abbreviations, the basic perturbation equation (37) is 
written as 

L + [2 + (1 - 13) gB] H5m - AnG.ffPmSm = , (44) 
with an effective gravitational constant 



Geff = G il + a 



1 - Qr 



(1-/3) 
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1+9 



5-2^UUA + fi1-4-^^ 



(45) 



Recall, that the interaction also explicitly enters the "friction" term in addition to its infiu- 



ence on the Hubble rate itself. This additional term in the factor that multiplies 6m in (44) 
vanishes only for (3 = 1. Changing to a as independent variable. 



Sm = S'mCtH , Sm = H'^ 



where 5L = 4?, we obtain the final equation 



1 H' 

<rf/ _|_ _|_ s:/ 

^m ' ^m ' TT m 

a n. 



(46) 



C + U{a)5'^ + V{a)5m = 



(47) 



with 



and 



where 



U{a) 



+ A{a) + {I - (3) g{a)B{a) 



A 



3 H' 



V{a) 



2a? G 



Vtm = {^mo + K [a^exp(-aVa2) - exp(-l/(T2)] } 



1^ 
a3 m 



(48) 



(49) 



(50) 
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/V,min 


a 


/3 


7.457 


"•"J^' -1.403 


1 205+2-504 



TABLE II: Best-fit values for the parameters a and /3, based on the growth-rate data in [2T] and 



With f2mo = 1 — Kexp(— 1/0-2) ^|^2 _ ^^^^ realizes that f2m ~ 1 both for high redshifts 
a ^ 1 and in the long-time limit a ^ 1. 



Equation (47) is the central equation for the Newtonian perturbation analysis. The free 
parameters are iJo? ^ aiid o as well as a and /3. In the non-interacting limit and for the 
one-component case p = Pm one recovers the perturbation equation 



A" -I- s' 



2a2^™ 







(51) 



for the Einstein-de Sitter universe. Moreover, (51) is also the limit of (47) both for a ^ 1 
and for a ^ 1, since all the interaction terms are vanishing under these conditions. In 
general, the interaction influences the perturbation dynamics through all the K terms in 



the coefficients U{a) and V{a) in (47) which also includes the modification of the Hubble 
rate due to the interaction on the background level. To test this model, we performed a 
Bayesian statistical analysis, based on the growth rate data collected in [21] and those of the 
WiggleZ survey |22]. The left panel of Fig. |4] shows the dependence of 5m for different values 
for a and (3. The thick solid (blue) curve corresponds to the best-fit values of Table |ll| 
visualized in Fig. |3| The right panel compares the best-fit curve with the corresponding 
behavior for an Einstein-de Sitter universe and with the ACDM model. The growth of the 
matter perturbations is reduced compared with its ACDM counterpart for values a ~ 1. For 
larger values of the scale factor the perturbations grow again, whereas they stay constant 
for the ACDM model. If the recent WiggleZ data are not included, there exists a maximum 
in 5m{a) close to the present time, similar to some of the curves in the left panel. 

The different growth of matter perturbations has been used in the literature primarily to 
discriminate between a GR-based behavior and alternative theories of gravity pTl l32lflT] . 
However, also the impact of interactions on the perturbation growth has been investigated 
[271 ED Il2] . It is convenient to introduce the growth rate 

d In 5m. 



f 



din a 



(52) 
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pap 



FIG. 3: Statistical analysis for the parameters a and /3, based on the grow-rate data in [21j and 
|22j . Left panel: contour plot for the best-fit values of Table [TH Center and right panels: probability 
distribution functions (PDFs) for a and /3, respectively. 




a 



FIG. 4: Fractional density perturbation as a function of the scale factor. Left panel: perturbation 
behavior for different values of a and /3. Right panel: comparison between our model with the best- 
fit parameters of Tables [I] and [n| and the ACDM model. The straight line shows the corresponding 
increase for the Einstein-de Sitter model. 



in terms of which the basic equation (47) takes the form 

df 



dhia 

The last equation can also be written as 



+ + [aU{a) -l]f + a' V{a) = . (53) 



+ f + [aUia)-l]f = l^n^. (54) 



din a 2 G 

In general, the effective gravitational constant Gg// differs from G due to the interaction 



terms. The first correction term (parameter a) in (45) describes the direct coupling to the 
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ff=-0.36,^=1.21 



a=-0.36,/i=1.21 : 



FIG. 5: The ratio G^ff/G as function of the scale factor for different values of a and /3. The solid 
blue lines correspond to the best fit. Left panel: evolution from the past until a = 1.2. Right 
panel: future evolution. Recall that the dark-energy density becomes negative for a > y^2/3a ~ 6, 
but is exponentially suppressed. The ratio G^ff/G approaches the asymptotic limit Gf.ff/G = 1 
after passing through a second maximum at a w 19. 
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FIG. 6: The quantity gB as function of the scale factor for the best-fit values of a and /3. Left panel: 
evolution from the past until a = 4. Right panel: future evolution. Recall that the dark-energy 
density becomes negative for a > y^2/3a f« 6, but is exponentially suppressed. The quantity gB 
approaches the asymptotic limit gB = after passing through a minimum at a ~ 17. 



dark-energy fluctuations, the second term encodes the modifications due to the perturbed 
interaction quantity Q. For /3 = 1 the second contribution vanishes: an interaction which 
is only operative in the background does not modify the effective gravitational "constant" . 
In such a case there is an influence of the interaction on / only through the quantity A{a) 
in the coefficient U{a) (cf. Eqs. (48) and (41)). The effective gravitational "constant" G^ff 



approaches G in the early matter-dominated phase a <^ 1 and also for a ^ 1, where the 
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matter dominates again. The behavior of Geff is shown in Fig. [5j The left panel shows the 
evolution from small values of a in the past until a = 1.2. The future behavior is depicted 
in the right panel. Recall that for a > a/2/3(t ~ 6 the dark energy density becomes 
negative but it is exponentially suppressed. After passing through a second maximum, in 
the far future a ^ 22 one has Geff = G again. Fig. |6] shows the corresponding behavior for 
the best-fit values of the quantity gB which modifies the friction term in the perturbation 
equation. In Fig. [T] the growth rate is contrasted with the observations summarized in [21] 
as well as with those of [2Zj and with the ACDM model. Around the present epoch (a ^ 1) 
the deviation from the Einstein-de Sitter value is larger than that for the ACDM model. 
This corresponds to the slower growth of 5m (o) for values of the order of a 1 in Fig. |4| 
For larger values of a, however, 5m (o) continues to grow while one has 5m (a.) = const for the 
ACDM model. 



f(z) 




FIG. 7: Dependence of the growth rate f{z) on the redshift z. The data are taken from |21j (blue 
data points) and [22J (red data points). The shaded region denotes the la level, indicating a large 
dispersion. 



16 



IV. RELATIVISTIC PERTURBATION THEORY 



A. General relations 

We assume that the cosmic medium as a whole can be described by the energy momentum 
tensor of a perfect fluid (neglecting anisotropic stresses and energy fluxes in the rest frame), 

Tik = pUiUk + phik , = 0) (55) 

where hik = Qik + UiUk and gikU^u^ = — 1- The quantity denotes the total four- velocity of 
the cosmic substratum. Latin indices run from to 3. 

We assume a split of Tik into a matter component (subindex m) and a dark energy 
component (subindex x), 

Both these contributions are assumed to have a perfect-fluid structure as well, i.e., 
and 

= p,u\ul + p.K' , h^^ = g'' + v'A ■ (58) 
Next, we admit an interaction between both components according to 

r^,k = Q\ K\k = -Q'- (59) 

Then, the separate energy-balance equations are (cf. [27]) 

- UmiT^ -k = Pm,aKn + {Pm + Pm) = -UmaQ"" , {gikK^u'^ = -1) (60) 

and 

- MxiTf = p^. a< + {p^ + p^) = u^aQ" , {gikv^x^l = -1) • (61) 

Each component has its own four-velocity. The quantities and G^; are defined as G^ = 
u^.^ and G^,. = respectively . For the homogeneous and isotropic background we assume 
ul, = u°;. = u"". Likewise, we have the momentum balances 



^mi^m ■,k ~ (Pm + Pm) O'm + Pm,ihm — KfniQ^ (^2) 
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and 

KirH",, = (p. + p.) < + p.,K^ = -K,Q\ (63) 

where aj^ = mJ^.^m^ and a" = tt^.^M^- The source term is spht into parts proportional 
and perpendicular to the total four-velocity according to 

Qi = u'Q + , (64) 
where Q = —UiQ"^ and = h\Q°-, with UiQ^ = 0. 



B. The case px = —px 

The contribution T!^ is supposed to describe some form of dark energy. In the simple 
case of an equation of state px = —px, where px is not necessarily constant, we have 

= -pj' ■ (65) 



In the background, the balances (60) and (61) take the forms 



p^ + ?,Hpm = Q° (66) 

and 

P. = -Q° , (67) 

respectively. Denoting first-order perturbations again by a hat symbol and recalling that for 
the background = = is valid, the perturbed time components of the four- velocities 
are 

Uo = u'^ = ut = ul = ^9oo ■ (68) 



According to the perfect-fluid structure of both the total energy-momentum tensor (55) 



and the energy-momentum tensors of the components in (58), and with = = 
in the background, we have first-order energy-density perturbations p = Pm + Px, pressure 
perturbations p = Pm + Px = Px and 

Ta = T^a + Txa =^ {p + P) Ua = PmUma + {Px + Px) Uxa ■ (69) 

Greek indices run from 1 to 3. For px = —px it follows 

Px = -Px P+P = Pm Uraa = • (70) 
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Since the component m is supposed to describe matter, it is clear from (69) that the per- 



turbed matter velocity Uma coincides with the total velocity perturbation Uq. With = u^' 



up to first order, the energy balance in (60) (correct up to first order) can be written as 



On the other hand, the total energy balance is 



= -0(p + p) 



For the difference it follows that 



P ~ Pm = {P - Pm) aU"' = ^aQ'^ 



(71) 



(72) 



(73) 



Since, at least up to linear order, p — Pm = px, equation (73) is equivalent (up to the first 
order) to 

Px = Px,aU'' = UaQ" ■ (74) 



In zeroth order we recover (67). The first-order equation is (cf. (68)) 



Px + PxU^ = {UaQ" 



(75) 



Notice that ( 75 ) results from a combination of the total energy conservation and the matter 



energy balance. It has to be consistent with the dark energy balance (61). At first order, 
the latter becomes 

Px + PxU^ = {UxaQ") ■ (76) 

This means that 

{uZQ^) = , (77) 

i.e., the projections of along u^a and along Ua coincide. Explicitly, 



(78) 



In a next step we consider the momentum balances. The total momentum conservation is 
described by 

htT% = (p„ + + p,) a« + = . (79) 

With Px = —px we have 

h'lT% = pma'' + = . (80) 
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Using u"^ = again, the momentum balance (62) for the matter component becomes 



(81) 



Notice that we have only used the total momentum conservation and the matter momentum 



balance. The momentum balance (63 ) of the dark energy degenerates for the case px = —px- 



It does not describe any dynamics. 

Again we introduce the fractional perturbation S„ 
energy balance takes the form 

Pm 



— in terms of which the first-order 

pm 



e 



Q 



With 



tta — Uafl — -g00,a — ^a,0 + 



and Ua = v^a one finds from the momentum balance (cf. [43j) 



V,a + 4>,a = - — [Px,a + PxV,a] 
Pm 



V + 



Pr, 



[Px + PxV] 



(82) 



(83) 



(84) 



At this stage the relation to the previous Newtonian treatment becomes evident. Neglecting 



the perturbations on the right-hand side of the momentum balance in (84) and replacing 



f — )■ Of we recover the result (36) of the Newtonian analysis. Neglecting the terms multiplied 



by in (82) and identifying G accordingly, the non-relativistic relation (35) is reproduced 



The dynamics of the expansion scalar is determined by the Raychaudhuri equation which 
in our case takes the form 



e + - al + A-kG (p + 3p) = . 



(85) 



In the background, H = —AnG {p + p) = —AnGpm is valid. The term in the Raychaud 
huri equation becomes at first order 

1 



a'^Pn 



{Ap^ + p^Av) 



(86) 



where A is the three-dimensional Laplacian. 

For the perturbed time derivative of the expansion scalar we have 



e = e + = e - e</. 



(87) 



20 



Consequently, the first two terms of the perturbed Raychaudhuri equation are 



0_ 90+ 06 . 



(88) 



For the derivative of the expansion scalar O = —127iGpm is valid. For the perturbations of 
the term AttG [p + 3p] we find 



AnG [p + 3p] = AtiG {pmSm + Px^x + 3pa;) 



(89) 



C. The Gauge-invariant perturbation equation 



It is convenient now to introduce gauge-invariant quantities to describe the perturbation 
dynamics by 

(90) 



rc s- I rr r I Px ^ < • 

Om = Om^ V , b^ = d^^ f , p^=p^+ p^V 



Pr. 



Px 



as well as 



Qc = Q^Qy ^ and Q^ = Q + Qv . (91) 
The superscript c stands for comoving. All the symbols have their physical meaning on 



comoving hypersurfaces v = 0. From (84) with = c^p^, where Cs is the sound speed in 
the rest-frame f = 0, it follows that 



,2 P=^ rc 



• I J, 2 rx r 

V + (p = — c, — d, 

P' 



(92) 



Equation (82) takes the form 



Q 



XC I ^ xc I Ac I 2Pm Px rc Ac 
Pm Pm Pm Pm 



Since 



we may rewrite Q'^ as 



Q = -Px 



Q" = -px^l - px 



l + C 



2 Ar 
' pm 



6' 



(93) 



(94) 



(95) 



The gauge-invariant perturbation of the interaction quantity Q is determined by the dark- 
energy density perturbations and their first derivative. It follows that 



Pm Pm J Pm 



Q_ 

Pm 



(96) 
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In terms of the gauge-invariant quantities, the first-order Raychaudhuri equation becomes 

(97) 



In a next step we differentiate ( 96 ) and use 



Px 

pm 



and 



H 



_ Px_ 
pm 

-A-nGpm 



3H 



Pm Px 



(98) 



(99) 



as well as (97). The result is 

Q 



px 



Pr. 



-^6^ - 

X 

Pm 



+ 



Q_ 

pm 



3H'c 



2„2 





\2H^ 


^m- 


Pm) 


Pm 








\ R.' 


Pm 


V Prr 


I / Pm _ 


\P_X_ 




Q 

H— 


, Pm 




Pm 



5^ + ATiGpJ^^ + 



a?p, 



2 - 3c: 1 



px 
Pm 



6' 



(100) 



Now we change to the scale factor as independent variable according to (46). Applying 

H 2a p ' ^ ' 

and 



6 = a^H^ 



6"+ 1 



3pm\ S' 



2 p J a 



(102) 



the previous equation (|100|) can be transformed into 
'3 



+ 



3p^ _Q_ 
2^ 2 p Hp„ 



"m 



+ 



_3 3p^ 1 
2 2 p \p 



+ 



X 

pm 

\Prr. 
3px 



(6 - 3cl) 



px 



3p^ 
2 p 



2 + 



px 



Q_ 

m 

Q 



Q 



Hp„ 



Pm J Hpn 



+ 



Q 



Hp„ 



+ 



k'^ px 



2 p 



^'a^^ Pr, 



2-3cUl + - 



^2 ■ 



+ 3c^(5-^^ 



3p^ 
2 p 



(103) 



Here we have used 



APx 

am^p„ 



! P, 

•am^Pr, 



(104) 
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We have expressed the perturbations Q'^ of the interaction term in terms of the dark-energy 



perturbation 5^. and its first derivative according to (95). In order to obtain a closed second 



order equation for 5^, we assume a proportionahty between the dark-energy perturbations 
and the dark-matter perturbation by introducing 



(105) 



The parameter e quantifies the relative magnitude of the perturbations of the dark energy. It 
corresponds to the parameter a of the non-relativistic theory. To avoid misunderstandings, 
we use a different symbol here. In many studies dark-energy perturbations are neglected 
from the outset. However, this is strictly justified only for a cosmological constant. Neglect- 
ing these perturbation may lead to unreliable conclusions concerning the interpretation of 
observational data EI]- For the model dealt with in [27|, dark-energy perturbations were 
shown to be negligible on scales that are relevant for structure formation, but may play a 
role on super-horizon scales. 



With (cf. (40)) 



1 



Q_ 

Pm 



+ A + B 



Q 



Hpr. 



4- 



9 



we find 

where 
F(a) = 



8Z + F{a)bt + G{a)b^^ = 



a ( 1 + 

pm 



2^2 



and 



G{a) 



1 + 

Pm 



3 
2 



Px 

P 



3 
2 



gB 



2+P^]gB-{Q-3cl) 



+ 



A + B) gB -4 



9 



a^l+g 



- + A + B-3ci 
2 



1 + 



+ 15c- 



Pm 

p,r k 



gB 

2 1 



9 



a^l + g 



. — — 

2> 'p^am^ 



with 



Px 
Pm 



1 + K{a^ exp (-a2/a2) - §^2 exp (-1/^2)) 



(106) 



(107) 



(108) 



(109) 



(110) 
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and 



P 



HI 



-A 



(111) 



H 

It is only for e 7^ that we have a scale dependence. The quantities A{a), -jp and B, are 



given by (|4T|, (|42]) and (|43]), respectively, with (cf. ([29|) 

exp (— a^/a^) .ft' exp (— a^/cr^) 



-is: exp(- 1/(72) 1 - 1(72/^ exp(-l/a2) 
where we have used that Vtmo = 1 — ^exp(— l/cr2) (|cr2 _ i)_ 



(112) 



For e = we recover 



F{a) 



and 



G{a) 



3 3p^ 
"2 2 p 



1 



2 p 



0) 



- + A + B]gB-A^--^ 

2 J (T^ 1 + (yf 



0) 



(113) 



(114) 



These coefficients coincide with the corresponding coefficients of the Newtonian analysis for 
a = /3 = 0. Since Q'^ was related to 5^ and its first derivative in (95), the limit a = /3 = of 



the Newtonian theory corresponds to e = within the present relativistic analysis. Recall 



that ^ 
p 



lA. 



Our interest is the matter power spectrum, defined by = \Sm,k\ , where Sm,k is the 
Fourier component of the density contrast 6^. In order to choose appropriate initial condi- 
tions, we use the circumstance that at early times, i.e. for small scale factors a ^ 1, the 



equation (107) has the asymptotic Einstein-de Sitter form 

3 



A" -J- — — A' 

2(2 



2a2 



6r, 



0, 



(a < 1) , 



(115) 



which coincides with the corresponding equation of the ACDM model at that period. This 
allows us to relate our model to the ACDM model at high redshift. We shall benefit from 
the fact that the matter power spectrum for the ACDM model is well fitted by the BBKS 
transfer function |31]. Integrating the ACDM model back from today to a distant past, 
say z = 10^, we obtain the shape of the transfer function at that moment. The spectrum 
determined in this way is then used as initial condition for our model. This procedure is 
similar to that described in more detail in references [l5l HH]. 

In Fig. [8] we display the power spectrum, based on the data of the 2dFGRS project, for 
different values of e for = 1. The scale dependence via the coefficient G{a) in eq. (107) 
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is sensitive to the product ec^. The thick (blue) curve (e = —0.000023) represents the best 
overall fit. On large scales, however, the curve with e = 0.001 shows the better performance. 
For = 1 only a very small factor e is compatible with the data. Otherwise, there appear 
non-observed oscillations in the matter power spectrum which resemble a similar behavior 
in (generalized) Chaplygin gases |l3l HTj. The left panel of Fig. [9] provides an alternative 
illustration of the matter power spectrum. The right panel confirms that there is no scale 
dependence for = 0. There is no dependence on e either. The solid blue line just repro- 
duces the BBKS transfer function. We verified that also for small non-vanishing values of 
Cg, i.e. = 0.1, the resulting best-fit values for e are of the order 10~''...10~^. Although a 
constant e corresponds to a very rough approximation, these results indicate that fluctua- 
tions of the dark-energy component are small indeed on scales that are relevant for galaxy 
formation. On the other hand, if they were zero exactly, there would be no scale dependence 
at all. As the left panel of Fig. |9] shows, even a very small value of e, although considerably 
larger than the best-fit value, influences the spectrum substantially on larger scales. This 
demonstrates an increasing role of the dark-energy perturbations with increasing scale. The 
point here is that the much larger number of data for the smallest scales has more weight in 
the statistical analysis than the fewer data on larger scales. Obviously, the same constant 
value of (an almost vanishing) e which gives a correct description on the smallest scales is 
not adequate on the larger scales of the sample (left panel of Fig. [9]). Conversely, a good 
performance on large scales is not compatible with the observations on intermediate and 
small scales. Consequently, for a more advanced analysis, a scale-dependent e should be 
used. 

D. Non-adiabatic perturbations 

Perturbations in an interacting two-component system are necessarily non-adiabatic. The 
non-adiabatic part of the pressure perturbations is p — |p. For our special case the crucial 
quantity is 

P - = (C^ + 1) Pi + ^ f ^ - ^) = (C^ + 1) + - [p'ir.P. - plPm] . (116) 
p ^ ' p \pm Px) P 
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FIG. 8: Matter power spectrum for Cg = 1 and different values of e. The thick soUd (blue) curve 
(e = —0.000023) represents the best overall fit. On large scales, however the curve with e = 0.001 
shows the better performance. Larger values of e result in (non-observed) oscillations. The data 
are taken from 1231 . 



Introducing here the balances 
Pi = ^^P'n.^ we find 



-Q and pm = —<iHpm + Q as well as 6^ = e6^ 



p-Pp = Npl 
P 



where 



N = 



Q _^ ^px_ 

3Hpm Pm 



Q 



SHpr, 



;ii7) 



;ii8) 



Even for e = 0, i.e. without fluctuations of the dark-energy component, the interaction term 



induces a non-adiabatic contribution to the to total pressure perturbation, given by (30) 
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FIG. 9: Matter power spectrum for = 1 (left panel) and = (right panel). While the total 
best-fit value for = 1 is e = —0.000023, it is obvious that on larger scales the dashed curve with 
e = 0.001 gives a better description. This corresponds to the expectation that the dark-energy 
perturbations are more relevant on the largest scales. For = the resulting curves do not 
depend on e. The data are taken from 



with g from (112) and ^ from (110). The quantity characterizes the non-adiabatic part 



of the pressure perturbations. Its dependence on the scale factor is visualized in Fig. [TO] For 
a > \J\(y ~ 6 the dark-energy density becomes negative, albeit exponentially suppressed. 

approaches its asymptotic value = for a ^ 22 after passing through a minimum at 
a ~ 17. The behavior of iV(a) is almost independent of and e. 



V. SUMMARY 



We have tested a phenomenological model of transient accelerated expansion in which an 
interaction in the dark sector is constitutive for the cosmological dynamics. The interaction 
has both to cancel a "bare" cosmological constant and, at the same time, to generate a 
phase of accelerated expansion. While the detailed structure of the interaction was chosen 
for mathematical convenience, it admits an analytic solution of the background dynamics, 
we think that it can be used to discuss general features of transient acceleration models. We 
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FIG. 10: Dependence of the non-adiabaticity factor N on the scale factor. The non-adiabaticity is 



but this part is suppressed exponentiahy. The asymptotic value = is approached at a ^ 22 
after has passed a minimum Nmin ~ —4 at a ~ 17. 

reconsidered the background dynamics of this model and performed a statistical analysis 
based on the SNIa data of the Constitution sample. The model predicts a future minimum 
of the deceleration parameter q which afterwards switches to positive values again (Fig. |2]). 
The dark-energy density becomes negative for a ^ 6 but this contribution is exponentially 
suppressed. The direction of the energy transfer is from dark energy to dark matter for 
a ^ 8. It is reversed for a ^ 8 but the entire interaction term becomes exponentially 
suppressed as well. 

The perturbation dynamics of the model was investigated both on the Newtonian and 
on the GR levels. Using a simple parametrization for the perturbed interaction term in 
the Newtonian setting and including perturbations of the dark-energy component within a 
simple ansatz, we carried out a statistical analysis, using the growth-rate data in [2Ij and 
[22] . This allowed us to quantify, although with a large dispersion, the role of the perturbed 
interaction term and the contribution of the dark-energy perturbations. For a ^ 2 the 
fractional matter perturbations deviate stronger from the Einstein-de Sitter behavior than 
for the ACDM model (Fig. lil). This corresponds to a larger difference for the growth rate 



negligible for a <^ 1 and for a ^ 1. For a > 




^0" ~ 6 the dark-energy density becomes negative 
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f{z) as well. For a ^ 2, however, the matter perturbations continue to grow while they 
approach a constant value for the ACDM model. The effective gravitational constant Ge// 
deviates considerably from G for a > 1 but the ratio Gef f/G approaches unity again in the 
far-future limit (Fig. [sj. 

The relativistic perturbation analysis for the interacting two-component system was per- 
formed in terms of gauge-invariant quantities with physical interpretation in the comoving 
gauge. Perturbations in such type of systems are intrinsically non-adiabatic. To decouple the 
relevant perturbation equation we assumed for simplicity that the fractional perturbations 
of the dark energy are proportional to the fractional matter perturbations. The statistical 
analysis, based on the data from the 2dFGRS project reveals that the factor of proportion- 
ality is very small. In other words, dark-energy perturbations are small on scales that are 
relevant for structure formation. However, our analysis also shows that the considered data 
range is not adequately described by a constant factor. For the smallest scales we have much 
more data than for larger scales. Consequently, the small-scale date have a higher statistical 
weight than the fewer data on larger scales. On the smallest scales, the dark-energy fluctua- 
tions are irrelevant indeed. On the other hand it is obvious that on larger scales the overall 
best-fit curve does not provide a good description of the observations. A considerably larger 
(but still small) value of the mentioned factor shows a much better performance (Fig. |9]). 
This indicates an increasing role of the dark-energy perturbations with increasing scale, a 
subject that deserves attention in future research. 
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